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Abstract: Detecting change features of climate variables in arid/semi-arid areas is essential for 
understanding related climate change patterns and the driving and evolution mechanism between climate 
and arid/semi-arid ecosystems. This paper takes Inner Mongolia of China, a unique arid/semi-arid 
ecosystem, as the study area. We first detected trend features of climate variables using the linear trend 
analysis method and then detected their trend-shift features using the breaks for additive seasonal and trend 
method based on the time-series of monthly precipitation and monthly mean temperature datasets from 
1962 to 2016. We analyzed the different change features of precipitation and temperature on a regional 
scale and in different ecological zones to discover the spatial heterogeneity of change features. The results 
showed that Inner Mongolia has become warmer-wetter during the past 54 years. The regional annual mean 
temperature increased 0.4°C per decade with a change rate of 56.2%. The regional annual precipitation 
increased 0.07 mm per decade with a slightly change rate of about 1.7%, but the trend was not statistically 
significant. The warmer trend was contributed by the same positive trend in each season, while the wetter 
trend was contributed by the negative trend of the summer precipitation and the positive trend of the other 
three seasons. The regional monthly precipitation series had a trend-shift pattern with a structural 
breakpoint in the year 1999, while the regional monthly mean temperature series showed an increasing trend 
without a periodical trend-shift. After the year 2000, the warmer-wetter trend of the climate in Inner 
Mongolia was accelerated. The late 20% century was a key period, because the acceleration of the wetter 
trend in some local zones (I and II) and the alleviation of the warmer trend in some local zones (VII, VII 
and IX) occurred simultaneously. Moreover, the change features had a strong spatial heterogeneity, the 
southeastern and southwestern of Inner Mongolia went through a warmer-drier trend compared with the 
other areas. The spatio-temporal heterogeneity of the climate change features is a necessary background 
for various types of research, such as regional climate change, the evolution of arid/semi-arid ecosystems, 
and the interaction mechanisms between climate and arid/semi-arid ecosystems based on earth-system 
models in Inner Mongolia. 
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1 Introduction 


Climate change is one of the most important issues in global change research (Mahoney et al., 
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2004). The fifth "Intergovernmental Panel on Climate Change" report claimed that the global mean 
surface temperature has increased since the late 19 century, with a warming of 0.72°C (0.49°C to 
0.89°C) over the period 1951-2012 (Hartmann et al., 2013). It has predicted by several studies that 
severe storms and heatwaves prolonged drought and raised sea levels in the near future (Charles, 
2012; Dufresne et al., 2013; Leng et al., 2015; Mora et al., 2017) could dramatically influence the 
habitats of all life on the planet (Mora et al., 2017). Climate variables are the key parameters in 
many models based on the earth-atmosphere system. Accurately analyzing the change features of 
climate variables is not only essential for monitoring and forecasting climate change but also 
important for the sustainable development of terrestrial ecosystems. 

Arid/semi-arid ecosystems around the globe are especially fragile and extremely sensitive to 
external interference, such as climate change (Maestre and Cortina, 2002; Chen et al., 2012; Han 
et al., 2015; Zhang et al., 2016). Recent evidence has shown that dry lands would expand all over 
the world (Huang et al., 2016) and suffer more from the influence of climate change compared with 
humid areas (Huang et al., 2017). Therefore, it is a pressing need to improve the understanding of 
the global climate change and the interaction mechanisms between climate and arid/semi-arid 
ecosystems (Hu et al., 2015; Huang et al., 2015a, b; Wang et al., 2017). 

In northern and western China, arid/semi-arid areas are widely distributed and cover one-third 
of the Chinese land area. Inner Mongolia, an especially sensitive region in China, has one of the 
most unique arid/semi-arid ecosystems (Peng et al., 2013; Liu et al., 2018). Because of global 
climate change and rapid economic growth during the past decades, the regional ecosystem has 
changed greatly (Xu et al., 2016). The Chinese government has launched several ecological 
restoration programs in this area since 2000 (Chen et al., 2017; Shao et al., 2017), such as the 
"Three-North Shelterbelt Project", the "Grain for Green Project" and the "Beijing-Tianjin Sand 
Source Control Project", which have greatly influenced the relationship between humans and 
nature. Therefore, the Inner Mongolia of China has become one of the most attractive areas for the 
studies about arid/semi-arid ecosystem (Chen et al., 2013; Hu et al., 2015; Huang et al., 2015b; 
Zhang et al., 2018), such as climate change analysis (Peng et al., 2013), climate-vegetation 
relationship analysis (Bhandari et al., 2015; Tong et al., 2017) and land cover change (Guo et al., 
2014). 

Precipitation and temperature are two key climate variables controlling arid/semi-arid vegetation 
growth (Bao et al., 2010; You et al., 2010; Peng et al., 2013; An, 2015). They are also important 
parameters for the calculation of various arid indices (Huang et al., 2016) and for the construction 
of earth-atmosphere models (Fisher et al., 2017). Therefore, it is crucial to detect the detailed and 
accurate change features of the time-series of precipitation and temperature datasets in Inner 
Mongolia. Previous studies mainly focused on the overall trend analysis of these two variables and 
contrary results were found in different studies, e.g., Hu et al. (2015) found that Inner Mongolia 
became warmer-drier from 1961 to 2012, while others demonstrated a warmer-wetter trend (Geng 
etal., 2016; Hu et al., 2017; Qin et al., 2018). Meanwhile, the trend-shift features in the time-series 
of precipitation and temperature in Inner Mongolia are interesting but have rarely been analyzed. 
Hu et al. (2015) found that the increasing trend in global temperature has slowed down since the 
late 1990s. It was unknown whether temperature change in Inner Mongolia reflected the same 
pattern. Hu et al. (2017) also found a nonlinear change in precipitation in Inner Mongolia through 
observing the time-series change trajectory. The detection of climate trend-shift features in Inner 
Mongolia could not only reveal the detailed change patterns from an apparently consistent trend 
line, but also elucidate the coupled and response relationship between climate and arid/semi-arid 
ecosystems. It has been confirmed that there is a nonlinear relationship between climate variables 
and vegetation cover in arid/semi-arid ecosystems (Wu et al., 2017; Harrison et al., 2018). The 
detection of detailed and accurate change features, especially the nonlinear change features within 
a time-series of precipitation and temperature, is useful for understanding the nonlinear 
relationships between climate variables and vegetation ecosystems in Inner Mongolia. 

This study aims to detect trend and trend-shift features in precipitation and temperature over the 
period 1962-2016 in Inner Mongolia, China. The features are closely related to vegetation growth 
in arid/semi-arid ecosystems. We detect the overall trend in time-series of annual precipitation and 
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annual mean temperature datasets, and then detect the decade-scale trend-shift features in monthly 
precipitation and monthly mean temperature datasets. We also analyze the change features in 
different seasons and ecological zones to further describe the temporal-spatial heterogeneity of the 
change features on a local scale. Our results could help to understand the local detailed change 
mechanisms of key climate variables, and then to provide basic background knowledge to research 
on the relationship between climate and arid/semi-arid ecosystems in the local study areas. 


2 Materials and methods 


2.1 Study area 


The Inner Mongolia Autonomous Region of China (37°24'-53°23'N, 97°12'-126°04’E), with an 
area of about 1.18x10° km’, was chosen as the study area (Fig. 1). The mean annual precipitation 
varies from 50 mm in the northwest to 500 mm in the northeast while the annual mean temperature 
varies from 0°C to 8°C. Humid, semi-arid and arid zones are found from east to west. Soil types 
are black soil, chernozem, chestnut soil and sandy soil. Vegetation comprises forest, steppe and 
desert steppe (Hu et al., 2015). 
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Fig. 1 Distributions of the 11 ecological zones (I-XI) and the 46 meteorological stations (1-46) 


2.2 Data 


In this study, monthly precipitation and monthly mean temperature datasets were acquired from the 
website of the China Meteorological Data Service Center (http://data.cma.cn/). The datasets contain 
records of 46 meteorological stations located in the study area (Fig. 1) for the period from January 
1962 to May 2016 and the quality of the datasets was tested, showing no erroneous data. The 
missing data at the beginning or end of the time-series were replaced by the average values of those 
months in all the other years, while the remaining missing data were filled using the linear 
interpolation method (Hu et al., 2015). The missing data, overall, comprised less than 0.08% of all 
datasets. In each meteorological station, there are 653 points within each monthly precipitation and 
monthly mean temperature series. All the 653 points were used in analyzing decade-scale trend- 
shift features based on the breaks for additive season and trend (BFAST) method. Monthly mean 
temperature series were reconstructed to annual mean temperature series with 54 points from 1962 
to 2015 by averaging the monthly values within each year. Monthly precipitation series were 
reconstructed to annual precipitation series with 54 points from 1962 to 2015 by the accumulated 
monthly values within each year. The annual mean series were used to analysis the overall trend 
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based on linear trend detection method. The data in 2016 were not involved in reconstructing the 
annual mean series because there were not enough samples in this year. The series of different 
seasons were reconstructed by accumulating the value of different months within each year, i.e., 
spring (March to May), summer (June to August), autumn (September to November) and winter 
(December to February in the next year). These seasonal accumulated time series temperature and 
precipitation datasets were used to analysis seasonal trend features. 

The natural environment of Inner Mongolia has a strong spatial heterogeneity. The spatial 
heterogeneity analysis was based on a 1:5,000,000 ecological zone map of Inner Mongolia. The 
study area is divided into 11 ecological zones (Table 1; Fig. 1) according to climate, soil, landform 
and vegetation conditions (Xu et al., 2016). As the distributions of meteorological stations in the 
zones of VII and VIII are sparse and uneven, we merged these two small neighboring zones into 
one. 


Table 1 Meteorological stations in different ecological zones 


Ecological Area Characteristic Meteorological 
zong Landscape Vegetation station code 
I Da Hinggan Mountains Grassland Quercus mongolica 1,2,3, o 1, 
; Stipa krylovii, 
H Hulun Buir Plateau ie Mongolian Plateau, Cleistogenes 4,5,8,9 
PP caespitosa 
Ill Songliao Plain Grassland Ulmus pumila 13, 17, 19, 23 
IV Xilingol Plateau Hunshandak sandy land, Ulmus pumila 14, 15, 16, 18, 20 
grassland 
V West Liaohe River Plain Aal plain, irrigated Crops: corn and wheat 24, 25, 26, 29, 30 
agricultural area 
21, 22, 27, 28, 33, 
VI Ulanqab Plateau Red sand, steppe desert Stipa glareosa 34, 35, 36, 37 
Yinshan Mountains, : ee : ; 
VI&VIHI Hetao-Tumòchuan Alluvial plain, irrigated Crops: wheat, rice and 38, 39, 40 
j agricultural area, steppe cereal 
Plain 

IX Ordos Plateau Steppe desert Cistus ladanifer 42,45 

x ba sad High Gobi, desert Ephedra 41, 43, 44, 46 
XI Lower reaches of Heihe Alluvial plain, oasis Populus diversifolia 32 

River Oasis forest 


2.3 Methodologies 


We first detected the trend feature of climate variables based on the time-series of annual 
precipitation and annual mean temperature datasets with 54 points, then detected the trend-shift 
features of climate variables based on the time-series with 653 points. The above detection 
processes were first conducted on a regional scale and then in different ecological zones to discover 
the spatial heterogeneity of the change features of precipitation and temperature. 

2.3.1 Detection of trend features 

The linear regression detection method (Hu et al., 2015) was used to determine the trend of each 
climate variable (a=0.05). The trend slope can be defined by the least-squares estimation method 
as in Equation 1: 
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where Slope is the trend slope of the time-series, with Slope>0 means a positive trend and Slope<0 
means a negative trend; n is the number of years (n=54); x; is the value of the independent variable 
(i.e., years); and y; is the value of the dependent variable (i.e., annual mean temperature series and 
annual precipitation series) in the i” year. 


The significance of Slope can be determined using the F-test as in Equations 2—4: 
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where U is the sum of squares due to regression; Q is the sum of squares due to error; x is the mean 
value of the independent variable; and yis the mean value of the dependent variable (i.e., annual 
precipitation series and annual mean temperature series). The change rate of dependent variable 
over the study period, which directly expresses the magnitude of the trend change, can be 
determined using Equations 5 and 6: 


ee (Slope x x, +a) —(Slope x x, +a) ee (5) 
Slope x x, +a í 
a= y-Slopexx, (6) 


where Rate is the rate of change; a is the intercept of regression equation; and xı and x, are the start 
and the end values of the independent variables during the time-series, respectively. 

2.3.2 Detection of trend-shift features 

The trend-shift feature is a nonlinear feature, which can be detected by time-series breakpoints 
analysis method. A trend-shift was determined only when statistical structural changes occurred in 
the time-series datasets. In this study, we used BFAST method to detect the trend-shift features in 
the monthly precipitation series and monthly mean temperature series based on the BFAST package 
in R software (R Core Team, 2017). BFAST is a multiple breakpoint detection method based on the 
trend-seasonal decomposition method and structural change tests, which is developed from 
econometric models (Zeileis et al., 2001, 2005; Zeileis, 2005). It has been successfully used to 
detect trend and seasonal structural changes in the normalized difference vegetation index time- 
series datasets (Verbesselt et al., 2010a, b) and has been proved to be a powerful method to detect 
structural changes in time-series models. The method has been proved less sensitive to cyclical 
changes in time-series (Verbesselt et al., 2010a; Jong et al., 2013). BFAST used iterative strategies 
to gradually find the breakpoints based on a trend-seasonal decomposition model and structural 
change detection model (Zeileis et al., 2001, 2005; Zeileis, 2005; Jong et al., 2013). The final output 
model was determined by a best-fit model, i.e., the Bayesian information criterion (BIC). The 
results of BFAST do not only reflect whether the trajectory has the trend-shift features or not, but 
also the confidence interval and mutation date if the trend-shift really does occur. In this study, the 
significance level was 95% (a=0.05). 

The BFAST method is different from traditional decomposition methods, such as seasonal-trend 
decomposition based on locally weighted regression (STL), which is only applicable to a time- 
series without continuous breaking signals (Cleveland and Cleveland, 1990). The BFAST method 
used an initial seasonal model resulting from STL to start the model as in Equation 7, and then 
replaced the piecewise trend and seasonal models through iteration. The iteration continues until 
the number and location of breakpoints does not change in the next iteration. 
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where y; is the time-series observations; a; is the intercept; a2 is the slope of the trend; 7; is the 
amplitude; 6; is the phase (i.e., season); k is the harmonic terms; f is the frequency (e.g., f=12 in the 
monthly precipitation series and monthly mean temperature series); and & is the unobservable error 
at time ¢, i.e., standard deviation (ø). 

The piecewise trend or seasonal model in BFAST was formed using the structural change 
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detection model after Equation 7 and was rewritten as Equations 8—10: 


Y. =X; B+E, (8) 
x, = {1,1,sin(2alt/f),cos(2zlz/f),...,sin(2nkt/ f ),cos(2akt/ f )}" P (9) 
B={a,,a,,7,C0s6,, 7sind,,...,7,c0S6,,7,8ind, Y, (10) 


where x; and £ are the independent variable and the regression coefficient in the standard linear 
regression model, respectively. 

The moving sums (MOSUM) of the ordinary least squares (OLS) residuals process (OLS- 
MOSUM), a generalized fluctuation test (Zeileis et al., 2001, 2005; Zeileis, 2005), was used to 
determine the breakpoints of structural change. The regression parameter / can be estimated and 
tested by the least-squares estimation. The OLS-MOSUM process is the sum of a fixed number of 
residuals in a temporal window whose size is determined by the bandwidth parameter. The window 
moves over the whole time-series and the OLS-MOSUM process is defined as Equation 11: 


where MO, is the OLS-MOSUM process; B is the estimated parameter f} from a stable historical 


period @1, ..., n; and A is the bandwidth of MOSUM (A € (0, 1)). The parameter h in BFAST was 
set to 0.18, which indicated that we limited the minimum interval between two breakpoint years to 
about 120 months, i.e., 10 years. This predetermined value of A corresponded to several climate 
change studies which concerned changes on a decade scale. 

The hypothesis of structural change detection is that the value of MO; should be close to 0 and 
only fluctuates randomly during the window's movement if there is no trend-shift in the time-series. 
Once the trend-shift appears when the window is moving, the value of MO, will systematically 
deviate from 0. The trend-shift is determined when the absolute value of MO, exceeds a boundary 
that is asymptotically only crossed with a 5% probability under structural stability (Verbesselt et 
al., 2012). 

The BIC was used to select the best-fit model once all the possible models with different 
breakpoints had been set up in the previous two models. The BIC is calculated as in Equation 12: 


BIC =-2In(L)+k, (12) 


where L is the maximum likelihood of the model; and k is the number of parameters. We selected 
the model with the minimum value of BIC. The optimum number and position of the mutation 
points occur when the BIC minimum is satisfied. At this time, the overall equation-fitting residual 
and the model complexity achieve the best balance (Bai and Perron, 2003; Verbesselt et al., 2010a, b). 


3 Results 


3.1 Trend features in Inner Mongolia 


3.1.1 Trend features on a regional scale 


Figure 2 shows the trend features of the annual mean temperature series and annual precipitation 
series on a regional scale in the study area during 1962-2015. The annual mean temperature series 
increased significantly, with a change rate of 56.2% and an increase in magnitude of 0.4°C per 
decade (a=0.05). During the same period, the annual precipitation series increased slightly, with a 
change rate of 1.7% and an increase in magnitude of 0.07 mm per decade. However, the positive 
trend in the annual precipitation series was not statistically significant (a=0.05). The results 
indicated that, on the regional scale, although annual precipitation remains relatively low with 
insignificant increases, the study area becomes warmer-wetter with temperature continually 
increasing during 1962-2015. 
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Fig. 2 Trend feature of the annual precipitation series and annual mean temperature series on a regional scale 
(left-top) and in different ecological zones from 1962 to 2015 


The anomaly analysis further showed that the lasting positive trend in the annual mean 
temperature series (Fig. 3) was contributed by the negative anomalies in 1962—1992 and the 
positive anomalies in 1993-2015 (Fig. 3). The coldest year was in 1969, with an annual mean 
temperature of only 3.2°C, while the warmest year was in 2007, with the annual mean temperature 
of 6.6°C. The coefficient of variations in the annual mean temperature series was 21.0%, which 
was higher than that in the annual precipitation series (14.9%). During the same period, although 
the overall trend in the annual precipitation series changed slowly and slightly (1.7%), there was 
strong inter-annual variations, with a coefficient of variations up to 14.9%. The wettest year was in 
1998, with an annual precipitation of 32.3 mm and a percentage precipitation anomaly increasing 
to 50.9%, while the driest year was in 1965, with an annual precipitation of only 16.5 mm. 
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Fig.3 Anomaly analysis of the annual precipitation series (a) and annual mean temperature series (b) on a regional 
scale from 1962 to 2015. The gray dotted lines in each bar chart denote a deviation from the mean value for the 
period. 


3.1.2 Trend features in different ecological zones 


Although there was a warmer-wetter trend feature on a regional scale, the results were varied when 
we analyzed the trend features of the annual precipitation series and annual mean temperature series 
in different ecological zones. Not all the ecological zones followed the same warmer-wetter trend 
as the overall regional feature (Fig. 2; Table 2). 

Trend features of the annual mean temperature series in each ecological zone and meteorological 
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station showed a statistically significant warmer trend. But the change rates of the annual mean 
temperature series in different ecological zones and meteorological stations had a strong spatial 
heterogeneity. Specifically, the meteorological stations 39 in Yinshan Mountains in the zone of VII, 
38 on Hetao-Tumochuan Plain in VIII, and 42 in the Ordos Plateau in the zone of IX, were the top 
three stations with the highest increasing rates in temperature with the increases of 0.56°C, 0.51°C 
and 0.52°C per decade, respectively. Meanwhile, meteorological stations 29, 30 and 26 in West 
Liaohe Plain in the zone of V had the lowest increasing rates, with the increases of 0.15°C, 0.22°C 
and 0.26°C per decade, respectively. Therefore, temperature increase was mainly concentrated in 
the southwestern (VII, VIII and IX) and southeastern (V) Inner Mongolia. 

Trend change of the annual precipitation series in different ecological zones and meteorological 
stations also had an obvious spatial heterogeneity. Moreover, some of them even showed an 
absolute opposite change trend. The zones of I, II, VI, VII, VIII and X had 29 meteorological 
stations distributed in the northeastern and northwestern Inner Mongolia. The trend features of the 
annual precipitation series showed the same wetter trend as that at the regional scale. 
Meteorological stations 2, 10 and 6 in the zone of I were the top three stations, with the highest 
increasing rates in annual precipitation, with the increases of 4.93, 1.17 and 1.14 mm per decade, 
respectively. Meanwhile, the zones of III, IV, V, IX and XI, with 17 meteorological stations, showed 
drier trends during 1962-2015. These regions are mainly located in the central and southern Inner 
Mongolia, especially in the southeastern and southwestern areas. Meteorological stations 26, 29 
and 25 in West Liaohe River Plain in the zone of V had the highest decreasing rates in annual 
precipitation, with the decreases of 1.20, 0.95 and 0.70 mm per decade, respectively. 

Trend analyses of the annual mean temperature series and annual precipitation series on an 
ecological zone scale indicated that 6 of the 11 ecological zones had the warmer-wetter trend as at 
the regional scale. The other 5 ecological zones went through a warmer-drier trend, within which 
the zones of V and IX were the typical hotspot areas. 


Table 2 Trend features of the annual precipitation series and annual mean temperature series in different 
ecological zones 


Ecological Changes of annual precipitation series per decade Changes of annual mean temperature series per decade 
aoe (%) (mm) (%) CC) 
I 11.38 0.76 153.62 0.36 
II 2.37 0.11 177.35 0.36 
Il —0.46 —0.03 46.58 0.39 
IV 0.41 —0.02 406.10 0.41 
V -11.21 —0.66 23.66 0.27 
VI 1.40 0.05 72.62 0.42 
VIL, VII 6.66 0.28 45.17 0.50 
IX —0.47 —0.02 44.15 0.45 
X 19.09 0.36 26.67 0.37 
XI —18.34 —0.10 34.17 0.49 
heat 1.72 0.07 56.16 0.40 


3.1.3 Seasonal trend features 


Seasonal trend features were also analyzed to discover more details in the changes of temperature 
series and precipitation series from 1962 to 2015 (Fig. 4). Figures 4e-f demonstrated that the 
positive trend in the regional seasonal accumulated temperature series was contributed by the 
statistically significant positive trend in every season. Temperature in the vegetation's non-growing 
season (October to April) was found to increase more dramatically than in the growing season (May 
to September), in the order of winter>spring>autumn>summer. 
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Fig. 4 Seasonal trend features of regional accumulated precipitation series (a-d) and seasonal accumulated 
temperature series (e—f) series in different seasons during 1962-2015 


The seasonal accumulated precipitation series showed obvious positive trends in spring, autumn 
and winter (autumn>spring>winter) (Figs. 4a, c and d). However, the seasonal accumulated 
precipitation series in summer decreased slightly (Fig. 4b). In Inner Mongolia, summer 
precipitation accounts for 60%—70% of the annual amount according to our datasets and a previous 
study (Zhao et al., 2012). The increased precipitation in autumn, spring and winter was offset by 
the decreased amount in summer, so that the annual increase was non-significant. Seasonal trend 
feature analysis in precipitation series and temperature series showed that the study area went 
through a warmer-drier summer and warmer-wetter spring, autumn and winter during 1962-2015. 


3.2 Trend-shift features in Inner Mongolia 


3.2.1 Trend-shift features on a regional scale 

Results of trend-shift detection in the monthly mean temperature series and monthly precipitation 
series on a regional scale are shown in Figure 5. The monthly mean temperature series showed an 
increasing trend without a statistically significant periodical trend-shift, while the monthly 
precipitation series had an obvious trend-shift pattern. The structural breakpoint was identified 
being in April 1999 (i.e., Time of BP(s) 448; Fig. 5) in the precipitation series, which separated the 
whole series into two sub-periods (a@=0.05). The trend change rates in the first and second sub- 
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periods were about 22.5% and 41.6%, respectively. The trend in the second sub-period were 
dramatically higher than in the first sub-period followed by a statistically significant drop in 
precipitation in the year 1999. This indicates that the warmer-wetter climate trend in Inner 
Mongolia was accelerated by the significantly increased trend-shift of the monthly precipitation 
series after the year 2000. 
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Fig. 5 Trend-shift features in the monthly precipitation series (a) and monthly mean temperature series (b) from 
1962 to 2016. Green lines denote piecewise trends during the study period, the red line denotes the magnitude of 
the structural breakpoints, and the blue line represents the confidence interval of the detected breakpoints. Time of 
BP(s) 448 means the 448" point in the time series is detected to be the break point. 


3.2.2 Trend-shift features in different ecological zones 


Although the regional monthly mean temperature series did not have a statistically significant 
trend-shift feature, there are trend-shift patterns in monthly mean temperature from different 
ecological zones. Table 3 shows that the zones of VII, VIII and IX had significant trend-shift 
features in 1996. It seems that the local trend-shift signals in these three zones were weakened or 
concealed after averaging the change in the other regions. The three neighboring zones were 
concentrated in the southwestern Inner Mongolia. However, the trend-shift pattern in Figure 6 
further indicates that the positive trend of the monthly mean temperature series slowed down after 
a temporary positive jump in 1996. 

Although the regional monthly precipitation series indicated a significant trend-shift in 1999, not 
all the ecological zones had the same change signals. It showed that only the monthly precipitation 
series in the zones of I and II went through the same significant trend-shift pattern in 1999 (Table 3). 
Therefore, zone I and its neighboring zone II, located in northern Inner Mongolia, were typical trend- 
shift areas because all the other areas of Inner Mongolia indicated no significant trend-shift. The 
same shift signals of the monthly precipitation series in the zones of I and II and in the whole study 
area indicated that the trend-shift in the whole study area was likely induced by the trend-shift that 
happened in northern Inner Mongolia during the past 54 years. Besides, it was noticeable that 
although zone I had a high increase in trend (Table 2), the positive trend of the monthly precipitation 
series in region I significantly increased since 1999 after a temporary negative jump (Fig. 6). 
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Fig. 6 Trend-shift of monthly precipitation series in zone I (a) and monthly mean temperature series in the zones 
of VII and VIII (b). Yt, the original time-series datasets; St, the seasonal component; Tt, the trend component; and 
et, the residual component. 


Table 3  Trend-shift features in different ecological zones 
Ecological zone Monthly precipitation series Monthly mean temperature series 

I 1999t-T - 

Il 1999t-t - 

Il - — 

IV - - 

y = 7 

VI - - 
VI, VII = 19961+1 

IX - 1996T+t 

xX = = 


Notes: — means that there is no significant trend-shift feature; 1999 and 1996 refer to the detected breakpoint years; f and } mean 
positive and negative trends ofa sub-period segmented by the breakpoint, respectively; — and + mean negative and positive jump directions 
in the breakpoint years, respectively. 


Therefore, the period 1996-1999 was a key period because the accelerated increase of the 
monthly precipitation series and alleviated increase of the monthly mean temperature series 
occurred simultaneously after the key period (Fig. 6). The results showed that the temperature and 
precipitation in these key years controlled or contributed to the spatio-temporal change pattern in 
local zones (1, H, VII, VIII and IX). 


4 Discussion 


4.1 Change features in Inner Mongolia 


Regional trend change features indicated that there was a warmer-wetter trend the study area during 
1962-2016. The positive trend in the annual mean temperature series was contributed to by the 
same positive trend found in each season (winter>spring>autumn>summer), while the positive 
trend in the annual precipitation series was contributed by the negative trend of the summer 
accumulated precipitation series and the positive trend of the other three seasons 
(autumn>spring>winter). After the year 2000, the warmer-wetter trend of climate in Inner 
Mongolia was accelerated. The results reveal unique climate change characteristics in Inner 
Mongolia compared with the other areas. 

The regional annual mean temperature series had a significantly positive trend, with an increase 
in magnitude of 0.40°C per decade, which was close to the results of Hu et al. (2015). However, it 
was obviously higher than the global warming rate of 0.14°C per decade and the warming rate of 
0.37°C per decade in the Tibetan Plateau (Yang et al., 2014), which is another unique arid area in 
China. This indicated that the arid/semi-arid area experienced a more accelerated warming trend 
than other areas around the world, with trend features possessing strong spatial heterogeneity in the 
world's different arid/semi-arid areas. The conclusion that temperature in winter increased faster 
than in summer was consistent with previous studies (Borjigin, 2015), and February was found to 
have the strongest temperature increasing trend in a previous study (Hu et al., 2015). Although we 
did not find a trend-shift change in the regional monthly mean temperature series, we did find it in 
the local monthly mean temperature series, in 1996 around Ordos Plateau (zone IX), Yinshan 
Mountains (zone VII) and Hetao-Tumochuan Plain (zone VIII). The change pattern in these areas 
corresponded to the previous conclusions that the increasing trend in global temperature has slowed 
down since the late 1990s (Hu et al., 2017). 

With regard to the change in the regional annual precipitation series, although the overall trend 
in the annual precipitation series changed slowly and slightly with a change rate of 1.7%, there was 
a strong inter-annual variation, with a coefficient of variation up to 14.9%. Because the 
precipitation series had a strong inter-annual variation, with a relatively continuous dry period from 
2000 to 2010, several studies based on precipitation data prior to 2010 resulted in a decreasing 
trend in the precipitation series (Hu et al., 2015), which was also proven by our analysis. The strong 
inter-annual variation in the annual precipitation series might be the reason for contrary results in 
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previous studies based on different study periods. The negative trend in the summer accumulated 
precipitation series was also confirmed by other studies in Inner Mongolia (Geng et al., 2016; Hu 
et al., 2017), and this feature was consistent across the whole Mongolian Plateau in northeastern 
Asia (Qin et al., 2018). We further analyzed the regional summer accumulated precipitation series 
and the results showed no trend-shift features. The detected trend-shift features in the regional 
precipitation series were consistent with previous results indicated by observing the time-series 
change trajectory (Hu et al., 2017). 

The spatio-temporal heterogeneity of the temperature series and precipitation series begs more 
careful attention. After the late 20" century, the warmer-wetter trend of climate in Inner Mongolia 
accelerated. However, the West Liaohe River Plain (zone V) and Ordos Plateau (zone IX), located 
in the southeastern and southwestern Inner Mongolia, respectively, went through a warmer-drier 
trend different from that in other zones. More attention should be paid to water management for 
sustainable development of the environment in areas with more fragile ecosystems, such as the 
Horqin Sandy Land in the north part of zone V and the Mu Us Sandy Land in zone IX. In other 
pastoral areas, such as the Xilingol Plateau, although the precipitation decreased slightly, the 
obvious and dramatic temperature increase might have accelerated regional evaporation and 
somehow influenced the typical, widely distributed, warm-temperate steppe. In some traditional 
agricultural areas, such as the southern zone V, the decreased precipitation may have decreased the 
precipitation-related corn yield. 


4.2 Importance of capturing nonlinear climate change 


The nonlinear change pattern of climate variables on regional and ecological zone scales was an 
important piece of data for many scientific studies focusing on arid/semi-arid areas. The trend 
change feature can only indicate an overall change direction of the climate variables, while 
capturing nonlinear climate change can profile more details hidden in an apparently consistent trend 
line. 

The arid index (Huang et al., 2016) and the standardized precipitation evapotranspiration index 
(Yao et al., 2018) can describe and differentiate the evolution of arid and semi-arid ecosystems. 
Precipitation and temperature are the key parameters, and nonlinear trend-shift features in the 
monthly precipitation series and monthly mean temperature series may influence the values 
calculated using these indices in some specific years, or even the overall trend. Besides, the coupled 
and response relationship between climate drivers and arid/semi-arid ecosystems may change 
because of the trend-shift features hiding in the monthly mean temperature series and monthly 
precipitation series. Our results indicated that the regional precipitation-vegetation relationship 
may change after the structural break year of the precipitation series. Trend-shift features in 
precipitation may provide a new perspective for studies on arid ecological resilience, which focus 
on the nonlinear dynamics observed in arid ecosystems (Holling, 1973; Wu et al., 2017). 

Trend-shift change features in the monthly precipitation series and monthly mean temperature 
series are also important for earth-atmosphere system models. The complicated interaction between 
climate, human activity and vegetation is critical for almost all studies based on earth-system 
models. Nonlinear changes in climate factors may influence the reliability of some models because 
they assume the stability of the climate-vegetation relationship. For example, in the residual trends 
model (Li et al., 2012; Xu et al., 2018), which is designed for discriminating between different 
drivers of land cover change, ignoring the trend-shift of time-series, climate variables may result 
in an overestimation or underestimation of the contribution of climate when discriminating between 
human-induced and climate-induced vegetation changes. 


4.3 Limitations and further research 


This paper detected trend and trend-shift change features in the temperature series and precipitation 
series without analyzing the reasons for the changes. Why these changes happened, however, is 
also an important and interesting scientific problem. We believe that many factors may together 
influence the change results, such as the changes in other climate variables, interaction between 
regional temperature and regional precipitation, interaction between regional human activity and 


arid/semi-arid environment and atmospheric circulation on a continental or global scale. Further 
studies could focus on this fields when there are more related datasets available. 

We only used publicly available climate datasets from 46 meteorological stations. We believe 
that data from more dense meteorological stations may add further detail to our results. To do more 
specific research, further study could use more meteorological stations, or some publicly available 
gridded climate datasets, such as version 4 of the moderate resolution imaging spectroradiometer 
land surface temperature dataset, or gridded climate datasets released by Yang et al. (2014). 

As discussed in Section 4.2, the nonlinear interaction between climate factors and arid/semi-arid 
vegetation is an important topic. This paper did detect the nonlinear changes in climate factors but 
it did not further analyze how they influenced the change in vegetation and whether there was a 
nonlinear relationship between climate factors and arid/semi-arid vegetation. This was because the 
temporal accumulation of the regional vegetation dataset, i.e., remote sensing vegetation index 
(NDVI) or enhanced vegetation index, was not long enough to match the temporal accumulation 
of meteorological data in this study. However, the results from this paper provide important 
background information for further studies on the nonlinear relationship between climate and 
vegetation in arid/semi-arid areas. We believe that multivariate time-series nonlinear models and 
nonlinear system stability analysis methods, such as phase space analysis, may have the potential 
to aid further study of this issue. 


5 Conclusions 


This study took Inner Mongolia of China as a study area, initially using a simple linear trend 
detection method to detect trend features in the annual mean temperature series and annual 
precipitation series and then the BFAST method to detect trend-shift features in those series. The 
results were as follows: (1) A warmer-wetter trend was shown during the past 54 years from 1962 
to 2015. The regional annual mean temperature series increased 0.4°C per decade with a change 
rate of 56.2%. The regional annual precipitation series slightly increased 0.07 mm per decade with 
change rate of about 1.7%, but the trend was not statistically significant. (2) The wammer trend 
was attributed to the same positive trend in each season (winter>spring>autumn>summer), while 
the wetter trend was contributed by the negative trend of the summer precipitation and the positive 
trend of the other three seasons (autumn>spring>winter). It showed the study area went through a 
warmer-drier summer, and at the same time a warmer-wetter spring, autumn and winter. (3) The 
monthly precipitation series had a trend-shift pattern with a structural breakpoint in 1999 separating 
the whole series into two sub-periods. The trend change in the second period was dramatically 
higher than in the first period followed by an obvious drop in precipitation in 1999. The monthly 
mean temperature series showed an increased trend without periodical trend changes. Therefore, 
the warmer-wetter trend of climate in Inner Mongolia accelerated after the year 2000. (4) Local 
zones went through a warmer-drier trend compared with other areas, such as zone V and IX. The 
late 20" century was a key period because the acceleration of the wetter trend in zones I and II, and 
the alleviation of the warmer trend in zones VII, VIII and IX, occurred simultaneously. 

The result is important for analyzing arid/semi-arid climate change, elucidating the interaction 
mechanisms between climate and arid/semi-arid ecosystems, illuminating other arid-ecosystem 
processes and augmenting prediction studies based on earth-system models. Further studies could 
focus on the change reasons and use more meteorological stations or some publicly available 
gridded climate datasets to do more specific research. Multivariate time-series nonlinear models 
and nonlinear system stability analysis methods, such as phase space analysis, may be used to 
specify the nonlinear relationship between climate factors and arid/semi-arid vegetation. 
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